Bounded modes to the rescue of optical transmission 



Michael Sarrazin^'B and Jean Pol Vigr 



■0 



o 
o 

o 
O 
1^ 



o . 

"—I , 
> 

Oh- 

q : 
o 

•I— I ' 

Oh 



(N 
m 

d 
o 



X 



stood by the conservation of light momentum such as 



^ Laboratoire de Physique du Solide, Facultes Universitaires Notre-Dame de la Paix, 
61 rue de Bruxelles, B-5000 Namur, Belgium 

This paper presents a brief survey of the evolution of knowledge about diffraction gratings. After 
recalling some basic facts, historically and physically, we introduce the concept of Wood anomalies. 
Next, we present some recent works in order to introduce the role of bounded modes in transmission 
gratings. The consequences of these recent results are then introduced. 

This paper is a secondary publication, published in Europhysics News (EPN 38, 3 (2007) 27-31). In 
the present version, some additional notes have been added with related references. 
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In the beginning of the 19*'* century, T. Young and 
J. Fraunhofer built the first optical diffraction gratings 
and showed the role of optical diffraction in their be- 
haviour. Since that time, diffraction gratings have been 
used in a broad range of technological applications, from 
the earlier accurate optical spectrometers to the recent 
integrated optical devices. They are also used in a broad 
range of wavelength from X-ray to microwaves. In this 
context, diffraction gratings have been the cornerstones 
of many publications for almost two centuries. Surpris- 
ingly, nowadays the subject is far from being exhausted. 

Every student knows that a grating can spread a white 
light beam from an incandescent lamp into a continu- 
ous spectrum of colors. For a grating made with a one- 
dimensional lattice of thick-thin wires, this "rainbow" 
is duplicated many times, each one corresponding to a 
diffraction order m (See Fig.l for instance). These basic 
behaviours of diffraction gratings can be easily under- 



k, = k, + K 



(1) 

where kg is the momentum of the scattered light, k^ is 
the momentum of the incident light, and K is a vector of 
the reciprocal lattice of the grating such that K • a =27rm 
(m S Z), where a is the basis vector of the lattice. The 
general behaviour of optical grating results from this ba- 
sic assumption. For instance, in Fig.l, one considers a 
reflecting grating. An incident light i is then scattered 
into three diffraction orders. The zero order (m = 0) is 
known as the specular beam. In Fig.l, two other orders 
are considered, m — \ and m — —1. 
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FIG. 1: A reflecting optical grating, i is the incident light 
beam. The light beams labelled m are some of the diffracted 
orders and 9 is the angle of incidence. 
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I. WOOD ANOMALIES 

In 1902, R. W. Wood observed some unexpected pat- 
terns in the spectrum of light resolved by optical diffrac- 
tion gratings [1]. The Wood spectrum presented many 
unusual rapid variations of its intensity in certain narrow 
wavelengths bands. As a consequence, these effects un- 
explained by ordinary grating theory were named " Wood 
anomalies" . 

In 1907, Lord Rayleigh proposed an explanation of 
those anomalies [2]. Considering for instance the case 
of the Fig.l, one notes that the order becomes tangent 
to the surface grating before disappearing. In such a 
case, the light momentum component along the Oz axis 
of a diffracted order m, i.e. kz^m, becomes imaginary 
right after having been cancelled. Then, the diffracted 
order m = \ becomes evanescent (non-homogeneous or- 
der). When \<iz,m is real, the diffracted order is free to 
propagate along the Oz axis (homogeneous order). In 
the present case, as the order m = 1 turns to be non- 
homogeneous, its energy will be redistributed over the 
other orders. As a consequence, the diffracted beam in- 
tensity of the specular beam (m — 0), for instance, in- 
creases just as the diffracted order m ~ 1 vanishes. For 
an incidence 6 (see Fig.l), a diffracted order m becomes 
non-homogeneous for wavelengths greater than a specific 
value, the Rayleigh wavelength, for which an anomaly 
occurs. One then talks about a Rayleigh anomaly. 



Around 1938, U. Fano proposed another explanation 
[3], where the anomalies are related to a resonance ef- 
fect. Such a resonance conies from a coupling between an 
eigenmode of the grating and a non- homogeneous diffrac- 
tion order (one also talks about a resonant diffraction or- 
der), i.e. after it vanishes for wavelengths greater than 
the corresponding Rayleigh wavelength. 

From an experimental point of view, gratings can then 
show some very complex behaviour often hard to explain 
in details without the use of Maxwell's equations. In 
addition, solving Maxwell's equations in order to study 
the electromagnetic diffraction quickly becomes a diffi- 
cult problem, even with simple geometries. Except for 
some simple and powerful analytical works, such as those 
from Lord Rayleigh or U. Fano, it will be necessary to 
wait until the Sixties and the beginning of computer 
calculation to observe significant progress in the Wood 
anomalies comprehension. For instance, one can under- 
line the first fundamental numerical results of Hessel and 
Olincr in 1965 [4], and Maystre and Ncvicrc in 1977 [5,6]. 
Hessel and Oliner presented a wide study of the eigen- 
mode resonance role in Wood anomalies [4]. They have 
shown that depending on the type of periodic structure, 
the two kinds of anomalies i.e., Rayleigh anomalies or 
resonant anomalies may occur separately or are almost 
superimposed. 

Maystre and Neviere studied specific cases of such reso- 
nant anomalies [5,6]. For instance, they presented a wide 
study about "plasmon anomalies", which occurs when 
the surface plasmons of a metallic grating arc excited [5] . 
They also considered an anomaly that appears when a 
dielectric coating is deposited on a metallic grating, and 
corresponds to guided modes resonances in the dielectric 
layer [6]. 

It is important to underline that, in some recent pub- 
lications, some authors improperly reduce the Wood 
anomalies to the Rayleigh anomalies only. It is impor- 
tant to recall that the Wood anomaly concept covers the 
resonant anomalies (Fano anomalies) and the Rayleigh 
anomalies at the same time. In addition, the resonant 
part of the phenomenon can imply many kinds of reso- 
nances and not only the surface plasmons resonances. 

Amazingly, thirty years after the pioneer works, in the 
early 2V* century, the experimental and theoretical in- 
terest in these subjects has not dried up. 



II. OPTICAL TRANSMISSION IN A 
NANOSTRUCTURED SLAB 

In 1998, Ebbesen et al [7] reported the amazing prop- 
erties of the optical transmission of some optical gratings. 
These were made by drilling cylindrical holes in a thin 
metallic layer along a bidimcnsional lattice. The metal- 
lic layer was deposited on a glass substrate (see Fig. 2). 
These experiments renewed the motivation for investi- 
gating metallic gratings, in particular those with a two- 
dimensional lattice. The most attractive characteristic 
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FIG. 2: The typical bidimensional array of subwavelength 
holes. 

of their results was the peculiar wavelength dependence 
of the transmission (sec Fig. 3). The latter was defined 
as the ratio between the energy of the specular transmit- 
ted beam (m = 0) only and the energy of the incident 
beam. After these first experimental observations, the 
role of the thin metallic film surface plasmons (SPs) was 
put forward in order to explain the peculiar wavelength 
dependence of the transmission [7]. 

A surface plasmon on a plane surface is a non-radiative 
electromagnetic mode, associated with a collective exci- 
tation of electrons at the interface between a conduc- 
tor and an insulator. In this way, a surface plasmon 
cannot be excited directly by light, and it cannot decay 
spontaneously into photons. The non-radiative nature of 
SPs is due to the fact that interaction between light and 
SPs cannot simultaneously satisfy energy and momen- 
tum conservation laws. Moreover, the momentum con- 
servation requirement can be achieved by roughening or 
corrugating the metal surface, for instance. It is exactly 
the situation in those devices. In this case, the transmis- 
sion peaks were interpreted as SPs resonances [7]. Sur- 
prisingly, this approach seems to have been suggested 
independently of the knowledge of some earlier works 
about the refiecting one-dimensional gratings. Then, the 
first attempt to understand the Ebbesen experiments has 
missed some important earlier results of the optical grat- 
ings study. 

In this way, in these first interpretations, the exact 
role of SPs was not clearly assessed, and many questions 
remained to clarify completely the scattering processes 
involved in these experiments. 

Note that, some authors developed other explanations 
as an alternative to the SPs model. For instance, it was 
suggested that the transmission pattern could also be 
described as resonances of electromagnetic cavity modes 
(here the cavities are the holes in the studied gratings) 
[8]. More recently, T. Thio and H. Lezec also suggested 
a diffracted evanescent wave model [9]. 

In addition, some doubts against the SPs hypothe- 
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FIG. 3: Computed optical transmission of a thin gold film 
(200 nm) deposited on glass. 100 nm is set for the radius of 
each hole and the lattice parameter is taken to be a = 700 
nm. 
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FIG. 4: Diagrams of slab modes. The electrostatic potential 
V of modes is plotted as a function of the location in the 
slab. Above: antisymmetric and symmetric modes. Below: 
pure surface mode, with each interface decoupled from each 
other. 



sis relied on some experimental results about subwave- 
length hole arrays made in non- metallic films. Though 
SPs cannot exist, the transmission pattern was found to 
be very similar to that obtained with a metallic film [9] . 
Surprisingly, the very fact that the typical transmission 
pattern can be observed even in non-metallic systems 
convinced some authors to fully reject the SPs hypoth- 
esis and rather consider models involving nonresonant 
evanescent waves diffraction. To our knowledge, this 
point of view is not supported by recent results. In fact, 
the SPs must be replaced by other kinds of eigenmodes 
in the nonmetallic cases. 

As a consequence, despite alternative theoretical inter- 
pretations, a large experimental consensus, many theo- 
retical results, tend to prove that the SPs interpretation 
was correct in the metallic films cases, and must be ex- 
tended to a bounded mode interpretation in the most 
general cases [10-12]. It is only the transposition of the 
1960's results about one-dimensional reflection gratings 
to bidimensional transmission gratings. 

Let us clarify this. 



III. BOUNDED MODES AND TRANSMISSION 

Basically, the gratings considered here are simply 
drilled slabs. In this bounded mode, i.e. an 

eigenmode of the grating, is an electromagnetic field con- 
figuration such that the field is trapped along the slab. 
The mode can propagate along the slab, but not else- 
where. The amplitude of the field decays then, exponen- 
tially from each side of the slab. Typically, two amplitude 
patterns exist, which correspond to symmetric and anti- 
symmetric modes. As the thickness of the slab increases, 
in some cases, the modes can refer to pure surface modes 
(Fig. 4), i.e. both slab interfaces are decoupled. Surface 
plasmons are such modes for instance. 



As previously suggested, bounded modes cannot be ex- 
cited directly by incident light or decay spontaneously 
into photons. Moreover, the coupling between the out- 
side electromagnetic field and the bounded modes is al- 
lowed by the use of rough or corrugated surface, as in the 
present case. The eigenmodes, and their coupling with 
outside light, can be described mathematically as follow. 
As in quantum mechanics, a scattering problem can be 
treated in electromagnetism via the use of the scattering 
matrix formalism (S matrix). In such a representation, 
it is assumed that the electromagnetic field can be de- 
scribed by two super- vectors \Fin) and \Fscat) (with the 
use of "ket" representation), which correspond to the in- 
cident electromagnetic field that lit the studied device 
and to the field scattered by the device. For instance, 
each component of these vectors corresponds to a spe- 
cific diffraction order. Then, incident and scattered field 
are linked via the scattering matrix defined as 

\Fscat) = S{X) (2) 

so that S{X) contains all the physical information about 
the studied device, such as a diffraction grating. In ad- 
dition, it can be shown that S{X) can be explicitly de- 
fined from the Maxwell's equations. For instance, the 
simulations in our papers are based on such a method, 
which combines scattering matrix formalism with a plane 
wave representation of the fields. This technique provides 
a computation scheme for the amplitude and polariza- 
tion (s or p) of reflected and transmitted fields in any 
diffracted order. A fundamental key to the understand- 
ing of diffraction grating phenomenology is to define the 
eigenmodes of the studied device. Eigenmodes then obey 
the homogeneous equation 

S-\Xp)\Fe,ge„) = (3) 

in such a way that the eigenmodes can also be defined 
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FIG. 5: Computed optical transmission (black curve) and 
amplitude of the related resonant order (red curve) of a thin 
fictitious film (the permittivity is equal to — 25 + if ) deposited 
on glass. 500 nm is set for the radius of each hole and the 
lattice parameter is taken to be a = 1000 nm. A discontinuity 
appears at the Rayleigh wavelength Xuay 



by the poles \p of the scattering matrix. Then, the poles 
are the solutions of the equation 



det{5"i(Ap)} = 



(4) 



If we extract the singular part of >S'(A) corresponding to 
the eigenmodes of the structure, we can write S{X) in an 
analytical form as 



A — A,] 



ShW 



(5) 



This is a generalized Laurent series, where Rp are the 
residues associated with each pole Aj,. One notes that 
Xp = Xp +iXp are complex numbers, and that the imagi- 
nary part Ap can be linked to lifetimes of the eigenmodes. 
The divergent behaviour of the fractional terms under- 
lines the resonant character of the eigenmodes (one talks 
then about the resonant terms of the S matrix). Sh{X) 
is the holomorphic part of S{X) which corresponds to 
purely non-resonant processes. Ideally, in the vicinity of 
a specific pole, the holomorphic part can be assumed to 
be almost constant. Then, as the amplitude s(A) of a 
specific diffraction order is related to an element of the 
scattering matrix, the intensity / oc |s(A)| of an order 
can be easily derived. One then obtains 



If = 



Jn 



A 



/ 2 



{x-x^y + x^p 



(6) 



where A^ — X^ + iXl is the complex zero of s(A) in the 
vicinity of Ap in the complex plan. 

This function leads to a typical asymmetrical profile, 
a Fano profile, related to a Fano resonance. A typical 
Fano profile appears in Fig. 5 (black curve). Now, if the 
holomorphic part is considered as negligible, one gets the 
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FIG. 6: Diagrammatic representation of the processes respon- 
sible for the behaviour of the transmission properties. 
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{x-x^y + x^ 



(7) 



which gives a simple Lorentzian profile, related to a 
Breit-Wigner resonance (see for instance the red curve 
of Fig. 5). As shown in Fig. 5, for a given eigenmode (i.e. 
Ap imposed) , the value of A for which Ib b. is a maximum 
(i.e. A = A^) does not correspond to the location of If 
maximum or minimum (i.e. (a) or (b) in Fig. 5). In order 
to illustrate the physical mechanisms responsible for the 
behaviour described by the scattering matrix formalism, 
we represent the corresponding processes involved. 

In Fig. 6, circles A and B represent diffracting ele- 
ments e.g., holes. So an incident homogeneous wave i 
(in red) diffracts in A and generates a nonhomogeneous 
diffraction order (blue dashed line). Such an order is 
coupled with an eigenmode, which is characterized by a 
complex wavelength Ap. It then becomes possible to ex- 
cite this eigenmode, which leads to a feedback reaction 
on the nonhomogeneous order (one can talk then about 
a resonant order). This process is related to the resonant 
term of the scattering matrix. In Fig. 5, the red curve is 
the electromagnetic field amplitude calculated for such a 
non-homogeneous resonant order. The resonant diffrac- 
tion order diffracts then in B, and generates a contribu- 
tion to a homogenous diffraction order, in our case the 
specular beam. Thus one can ideally expect to observe 
a resonant profile, i.e. Lorentzian-like, for the amplitude 
Sr of the homogenous specular diffraction order that ap- 
pears in B. Nevertheless, it is also necessary to account 
for nonresonant diffraction processes related to the holo- 
morphic term ^^.(A). So the incident wave i, generates 
also a homogeneous order Sh (in green), which appears 
as another contribution to the specular beam. So, as one 
observes the specular transmitted beam, one observes the 
sum of two rates, i.e. Sr + Sh- So, the Fano profiles of 
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the transmission (but also of the reflection) result from a 
superposition of resonant and nonresonant contributions 
to the observed diffraction order [10-11]. As a result, in 
Fig. 5 one can observe the calculated transmission (black 
curve) in comparison with the non-homogeneous reso- 
nant order amplitude (red curve). 

IV. CONSEQUENCES AND OBSERVATIONS 

As shown in Fig. 3, the observed optical transmission 
through subwavelength hole arrays, exhibits a set of 
peaks and dips. Following the previously described ap- 
proach, we have shown in recent works that the trans- 
mission spectrum is better depicted as a scries of Fano 
profiles (see Fig. 3 and Fig. 5) [10,13,14]. These recogniz- 
able line shapes result from the interference of nonreso- 
nant transfers with resonant transfers, which involve the 
eigenmodes film, and evanescent diffraction orders. We 
could then point out that each transmission peak-dip pair 
is nothing other than a Fano profile. The appearance of 
an assymetric Fano line shape does not necessarily lo- 
cate the peak or the dip at the eigenmode resonance (as 
shown in Fig. 5), contrary to what has been suggested in 
some earlier works [7,8]. However, we have shown that 
the existence of the eigenmodes is a condition for the 
presence of the Fano line shapes. It can be noted that 
recent results by Genet et al confirm this description [11]. 
As a further outcome of this work it became clear that 
this kind of transmission spectrum, with its Fano line 
shapes, could also be obtained in a large context and that 
the generic concept of eigenmodes must substitute that 
of SPs only. These results implied two important ideas. 
First, according to the contributions of resonant and non- 
resonant processes in the Fano profile, eigenmodes can 
be associated with wavelengths closer to the peaks or 
to the dips. Second, it is possible to obtain transmission 
curves similar to those of metal films, by substituting SPs 
with guided modes or other types of polaritons. Many 
examples can be found, including highly refractive ma- 
terials defining guided modes or ionic crystals in the re- 
strahlen band defining phonon polaritons (the restrahlen 



band of an ionic crystal is a narrow domain of wavelength 
for which real part of the permittivity becomes negative 
like in metals). For instance, we have reported simula- 
tions of a device, which consists of arrays of subwave- 
length cylindrical holes in a tungsten layer deposited on 
a glass substrate [13]. Tungsten becomes dielectric on a 
restricted domain of wavelength in the range 240-920 nm, 
i.e., the real part of its permittivity becomes positive. So, 
whereas plasmons cannot exist, wc show that the trans- 
mission pattern is similar to that obtained in the case of 
a metallic film. Indeed, it was also shown that in this 
case, instead of SPs, guided modes are excited [13]. Nev- 
ertheless, these modes give less intense fields than those 
achieved with SPs. In this way the benefit for the trans- 
mission is less important than with SPs modes. In the 
same way, it was also demonstrated that the transmission 
profile of a chromium film, in the restricted wavelength 
domain (1112-1292 nm), where the dielectric constant is 
positive, should involve eigenmodes, which are not SPs 
or guided modes. In this case, SPs can be substituted by 
Brewster- Zennek (BZ) modes [14]. Lastly it should be 
noted that in both cases, these theoretical computations 
are supported by experimental results. 

As a consequence, it appears that bounded modes can 
act as a mediator of the optical transmission in bidimen- 
sional subwavelength hole arrays. Although it appears 
that the mechanisms involved in these devices have been 
implicitly known for a long time, it seemed important 
to us to recall it. Indeed, an important engine of the 
technological development is also the improvement of our 
theoretical knowledge of the phenomena with which we 
are confronted. 

Additional notes : 

- The theoretical results about chromium, film,s are cor- 
roborated by the experimental data of Ref. [15]. 

- It is interesting to underline that those works about sub- 
wavelength hole arrays have implications in many other 
dom,ains such as biosensor technologies or left-handed 
metam,aterials (see Ref. [16] for instance). 

- The reader can refer to papers in Ref. [17] for plas- 
monic physics and technology in the context introduced 
in the present paper. 
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